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ABSTRACT 

Galactic stellar halos have long been considered to contain fossil information on early dynamical 
and chemical evolution of galaxies. We propose that the surface brightness distributions of old stellar 
halos contain the influence of reionization on early formation histories of galaxies. By assuming 
that reionization significantly suppresses star formation in small subgalactic clumps virialized after 
reionization redshift (z rc i n), we first numerically investigate how structural and kinematical properties 
of stellar halos formed from merging of these subgalactic clumps depend on z rc j on . We then discuss 
what observable properties of galactic stellar halos offer us the fossil records of reionization influence on 
hierarchical formation of halos based on the current results of numerical simulations. We particularly 
suggest that both the half-light radius of stellar halos and the slope of their surface brightness profile 
contain useful information on when star formation in subgalactic clumps were significantly influenced 
by reionization. By using the simulated surface brightness distributions of galactic stellar halos for 
models with different z r oion, we also discuss how wide-field imaging studies of extragalactic halos will 
help us to elucidate the epoch of cosmic reionization. 

Subject headings: Galaxy: halo - galaxies:evolution - galaxies: stellar content 



1. INTRODUCTION 

Structural, kinematical, and chemical properties of 
old, metal-poor stellar halos of galaxies have long pro- 
vided valuable information on early formation histories 
of galaxies, in particular, galaxies in the Local Group 
(LG) (e.g., Freeman 1987; Majewski 1993; Helmi et al. 
1999; Morrison et al. 2000; Ibata et al. 2001). For exam- 
ple, the observed correlation between chemical and dy- 
namical properties of Galactic halo stars or lack thereof 
have enabled authors to discuss whether the Galaxy was 
formed by a monolithic collapse within a few dynamical 
timescales (Eggen, Lynden-Bell, & Sandage 1962; Norris 
& Ryan 1991; Chiba & Yoshii 1998). Numerical simula- 
tions of Galaxy formation based on the cold dark matter 
(CDM) model (e.g., Bekki & Chiba 2000, 2001) have 
demonstrated that basic physical processes involved in 
the formation of the metal-poor, old stellar halo in a disk 
galaxy are described by both dissipativc and dissipation- 
less merging of subgalactic clumps and their resultant 
tidal disruption in the course of gravitational contrac- 
tion of the galaxy at high redshift: Understanding of star 
formation and chemical evolution histories in subgalac- 
tic clumps (i.e., progenitors of low-mass dwarfs) are key 
for better understanding of the origin of galactic stellar 
halos. 

Recently, several authors have extensively investigated 
the possible physical effects of reionization on galaxy for- 
mation, in particular, on star formation histories and 
physical properties of old stellar populations in low-mass 
dwarfs (e.g., Bullock et al 2000; Gnedin 2000; Susa & 
Umemura 2004; Grebel & Gallagher 2004; Willman et 
al. 2004). Some of these studies suggested that suppres- 
sion of star formation resulting from reionization can be 



a very important physical process for better understand- 
ing the observability and the early star formation activ- 
ities in dwarfs. These previous results, combined with 
those on stellar halo formation via hierarchical merging 
of low-mass subgalactic clumps (Bekki & Chiba 2000; 
2001), therefore imply that the influence of reionization 
on the formation of dwarf galaxies may be imprinted in 
the physical properties of old stellar halos in bright galax- 
ies. 

Although many observational and theoretical studies 
in the young universe have already provided important 
constraints on the cosmic reionization history (e.g., Gunn 
& Peterson 1965; Benson et al. 2001; Fan et al. 2003; 
Santos 2004; Wyithe & Loeb 2004), it remains yet to be 
explored in the present-day universe as to (1) whether old 
stellar halos of galaxies contain the fossil records of reion- 
ization history and (2) what observations of stellar halos 
can provide what clues to the possible influences of reion- 
ization on galaxy formation processes. Physical proper- 
ties of metal-poor, old stellar halos in galaxies beyond the 
LG have been recently investigated by the Hubble Space 
Telescope (HST) and large ground-based telescopes (e.g., 
Harris et al. 2000; Gregg et al. 2004). Furthermore, re- 
cent numerical studies provided some specific predictions 
on star formation histories of dwarf galaxies under the 
influence of reionization (e.g., Susa Sz Umemura 2004). 
Thus it is quite timely to discuss the above questions 
based on the comparison between observations and nu- 
merical simulations addressing the possible influence of 
reionization on the properties of galactic stellar halos. 

The purpose of this Letter is to show that wide-field 
imaging of galactic halos by large ground-based tele- 
scopes (e.g., Suprime-Cam on Subaru) can provide vital 
clues to the possible influences of reionization on early 
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formation history of galaxies. We first present the re- 
sults of our numerical simulations and stress how struc- 
tural and kinematical properties of stellar halos, which 
are formed from merging of small subgalactic clumps, 
depend on the epoch of reionization (z ro ion)- We then 
discuss (1) the origin of diverse stellar halo properties of 
galaxies of the LG in terms of reionization effects on halo 
formation and (2) the importance of ongoing and future 
wide-field optical imaging of the halos in better under- 
standing galaxy formation. Possible reionization effects 
on structural properties of globular cluster systems in 
disk galaxies are discussed in Bckki (2005). 

2. THE MODEL 

We simulate the formation of a Milky Way sized galaxy 
halo in a ACDM Universe with Q = 0.3, A = 0.7, 
Ho = 70 km s _1 Mpc -1 , and as = 0.9, and thereby 
investigate merging/accretion histories of subhalos that 
can contain low mass dwarfs. The way to set up ini- 
tial conditions for the numerical simulations is essen- 
tially the same as that adopted by Katz & Gunn (1991) 
and Steinmetz & Miiller (1995). We consider an isolated 
homogeneous, rigidly rotating sphere, on which small- 
scale fluctuations according to a CDM power spectrum 
are superimposed. The initial total mass (M t ), radius, 
spin parameter (A), and the initial overdensity Si of the 
sphere in the standard model are 6.0 x lO n M0, 30 kpc, 
0.08, and 0.26, respectively. Initial conditions similar to 
those adopted in the present study are demonstrated to 
be plausible and realistic for the formation of the Galaxy 
(e.g., Steinmetz & Miiller 1995; Bekki & Chiba 2001). In 
addition to this standard models, we investigate models 
with different M t and Si. 

We start the collisionless simulation at z sta rt (=30) and 
follow it till z cn d (=1) to identify virialized subhalos with 
the densities larger than 170p c (z), where p c {z) is the crit- 
ical density of the universe, at a redshift z. This 170p c (z) 
corresponds to the mean mass denstiy of a collapsed and 
thus gravitationally bound object at z (e.g., Padmanab- 
han 1993). The minimum number of particles within a 
virialized subhalo (N m i n ) is set to be 32 corresponding 
to the mass resolution of 3.8 x 10 7 M@. This number 
of 32 is chosen so that we can find a virialized object 
at a given z in a robust manner. For each individual 
virialized subhalo with the virialized redshift of z v ; r , we 
estimate a radius (rb) within which 20 % of the total 
mass is included, and then the particles within rb are la- 
beled as "baryonic" particles. This procedure for defin- 
ing baryonic particles is based on the assumption that 
energy dissipation via radiative cooling allows baryon to 
fall into the deepest potential well of dark-matter halos. 
Such baryonic particles in a subhalo will be regarded as 
candidate "stellar" particles to form stellar halos in the 
later dynamical stage, if the subhalo is later destroyed 
and baryonic particles initially within the subhalo is dis- 
persed into the galactic halo region. Thus, the present 
dissipationless models track the formation of stellar halos 
via hierarchical merging of subhalos, although the models 
are not adequate to the study of star formation histories 
in subhalos (as was done in our previous studies, e.g., 
Bekki & Chiba 2001). 

Previous theoretical studies have demonstrated that 
ultraviolet background radiation in a rcionized universe 
can significantly reduce the total amount of cold HI and 



molecular gas that are observed to be indispensable for 
galactic active star formation (e.g., Susa & Umemura 
2004). In order to investigate this suppression effects of 
star formation on the final structural properties of the 
simulated stellar halos, we adopt the following idealized 
assumption: If a subhalo is virialized after the comple- 
tion of the reionization (z Te i on ), star formation is totally 
suppressed in such a subhalo. Then, hypothetical bary- 
onic particles in the subhalos with z v - K < z Ic i on will not 
be identified as stellar halo particles in the later stage, 
but those in the subhalos with z v i r > z ro i on will be re- 
garded as progenitors of visible stellar halos. Recent 
WMAP ( Wilkinson Microwave Anisotropy Probe) obser- 
vations have shown that plausible z ro i on ranges from 11 
to 30 (Spergel et al. 2003; Kogut et al. 2003) whereas 
quasar absorption line studies give the lower limit of 6.4 
for z ro ion (Fan et al. 2003). Guided by these observations, 
we investigate the models with z re ion = (no reioniza- 
tion), 6, 8, 10, 13, and 15. The adopted picture of single 
epoch of reionization might well be somewhat oversimpli- 
fied and less realistic (e.g., Gnedin 2004), however, this 
idealized model can help us to elucidate some essential 
ingredients of the reionization effects on stellar halo for- 
mation. 

All the calculations have been carried out on the 
GRAPE board (Sugimoto et al. 1990). Total number 
of particles used in our simulations is 508686 and the 
gravitational softening length is 0.38 kpc. We note that 
later accretion of satellites at z < 1.5 is minor in the final 
structures of the simulated halos at z = 0, so the calcu- 
lation is ended at z onc i = 1 to obtain the dynamically 
relaxed halo structures. We used the COSMICS (Cos- 
mological Initial Conditions and Microwave Anisotropy 
Codes) , which is a package of fortran programs for gen- 
erating Gaussian random initial conditions for nonlin- 
ear structure formation simulations (Bcrtschinger 1995). 
In order to derive 2D distributions of B— band surface 
brightness (/^b) from the simulated density distributions 
of stellar halos, we assume that the stellar mass-to-light- 
ratio of 4.0 that is a reasonable value for stars with ages 
of ~ 13 Gyr and mctallicities of [Fc/H] ~ —1.6 (Vazdckis 
et al. 1996). 

3. RESULT 

Figure 1 shows that the projected radial density pro- 
files of stellar halos depend strongly on z TO ion hi such a 
way that they are steeper for the models with higher 
Zreion- This result implies that radial density profiles of 
galactic stellar halos can contain fossil information on the 
reionization epoch. The physical reason for this depen- 
dence is described as follows. For the model with higher 
2 re ion (=15), only subgalactic clumps that are virialized 
before z re ion can have old stars (i.e., no suppression of 
star formation by reionization) within it. They merge 
with one another in the early stage of the galactic col- 
lapse to form the first, single massive stellar halo around 
z ~ 2 with a steep density profile resulting from vio- 
lent relaxation of the merging. After that, subgalactic 
chimps that are virialized later than z rc ; on and thus do 
not contain stellar components merge with and accreted 
by the first massive halo. These lately merging subgalac- 
tic clumps are highly likely to be destroyed in the outer 
part of the galaxy. Consequently, if these subgalactic 
clumps contain stars initially within them, they are dis- 
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persed into the outer halo region and thus significantly 
flatten the radial density profile of the stellar halo. For 
the model with high z ro i on (=15), such later accretion 
cannot flatten the radial density profile of the stellar halo 
so that it can show a steeper density profile. Structural 
evolutoin of dark-matter halos due to hierarchical clus- 
tering process has already been demonstrated by Zhao 
et al. (2003), and the same process is at work for the 
-Zrcion-dependence of stellar halos in the present study. 

Figure 2 shows how surface brightness (^b) distribu- 
tions of galactic stellar halos in the B band depend on 
Zrcion- It is clear from Figure 2 that the model with 
higher z rc j on (=15) shows a sharp edge in its projected 
light distribution of the stellar halo in comparison with 
the model with lower z rc i on (=6). The ratio of the half- 
mass radius of the stellar halo to that of the dark matter 
halo is 0.15 for z ro i on = 15 and 0.71 for z rc ion = 6, which 
suggests that 2; re i on can control the compactness of the 
stellar halo distribution of a galaxy. Figure 2 demon- 
strates that the outer stellar halo (10 < R < 20 kpc) is 
much fainter in the model with z ro i on = 15 (/iB ~ 28 — 30 
mag arcsec~ 2 ) than in that with 2; ro i on = 6 (pb ^25 — 27 
mag arcsec~ 2 ). 

Figure 2 furthermore shows that a few stellar sub- 
structures can be seen in the outer stellar halo only for 
the model with higher 2; ro i on (=15), which implies that 
the observability of stellar halo substructures (or inho- 
mogeneity in luminosity distribution) can also depend 
on z ro ; on . Irrespective of z Ic i on , the stellar halos have 
prolate-triaxial shapes that follow the shapes of the dark 
matter halos in the present dissipationless simulations. 
Therefore, the /iB distributions of stellar halos depend on 
from which direction we observed the halos. It is accord- 
ingly noted here that the differences in shapes of stellar 
halos seen in Figure 2 are not due to z ro i on - Thus the 
way of radially falling /iB of stellar halos (i.e., presence 
or non-presence of outer sharp edges in /xb), compactness 
of the halos, and the average value of /zb can be regarded 
as fossil records of z ro i on . 

Figure 3 shows that /iB distributions of galactic stel- 
lar halos depend on the initial mass densities of galax- 
ies (5i) in the sense that the galactic stellar halo in the 
model with larger <5; (=0.39), shows higher p^ as a whole. 
This is because total number of subgalactic clumps that 
can be virializcd well before z ro i on is larger in the model 
with larger <5j. This result implies that if reionization can 
strongly suppress star formation for subgalactic clumps 
virialized after z ro j on , pb distributions of stellar halos 
can be remarkably diverse even for galaxies with simi- 
lar masses (yet different 5;) owing to the ^-dependence 
of the reionization effects. 

Kinematical properties of stellar halos in the models 
with a given mass depend weakly on z ro j on . For example, 
the maximum rotational velocity V m of stellar halos is 
appreciably larger in the model with 2; rc i on = 6 (Vm ~ 20 
km s _1 ) than that with z ro i on = 15 (Vm ~ 10 km s _1 ). 
This result of more rotation in the model with smaller 
•Zroion is described as follows. For the model with z ro i on = 
6, subgalactic clumps which are virializcd at 5 < z < 15 
and later merge with the galaxy contain a larger amount 
of orbital angular momentum. Therefore, the outer halo 
stars formed from tidal destruction of these clumps can 
have a larger amount of intrinsic angular momentum af- 
ter conversion of orbital angular momentum into intrinsic 



one in the merging processes. 

4. DISCUSSIONS AND CONCLUSIONS 

Previous observations suggested that the radial den- 
sity profile (p(r)) of the Galactic stellar halo can be 
well described as p(r) oc r -3 5 corresponding roughly to 
S(i?) oc i?~ 2 ' 5 in the projected profiles (e.g., Sommer- 
Larsen & Zhen 1990; Chiba & Beers 2000). Considering 
the present result that p(r) of galactic stellar halos in the 
models with higher z ro i on (=15) are as steep as r~ 3 - 5 and 
thus consistent with the observations, we suggest that 
the origin of the observed steep p(r) could be closely 
associated with reionization effects on the formation of 
subgalactic clumps of the Galaxy at a relatively higher 
rcdshift (~ 15). Furthermore, the simulated halos with 
sharply truncated outer distributions in the models with 
higher z ro i on suggest that the observed apparently sharp 
edge in the M31 stellar halo (e.g., Ferguson et al. 2002) 
might well result from such reionization effects. 

One of the most remarkable results of previous obser- 
vations for the M31 stellar halo is that it is dominated 
by moderately high-metallicity ([Fe/H] ~ —0.5) popu- 
lations and shows a high stellar density compared with 
the Galactic stellar halo (e.g., Durrell, Harris, & Pritchct 
2001; Reitzel & Guhathakurta 2001). The origin of the 
high metallicity, high-density stellar halo of M31 was dis- 
cussed in terms of the formation of M31's bulge via past 
major merger event (Bekki et al. 2003). We here sug- 
gest an alternative scenario that the observed differences 
in metallicity distributions and structures of stellar ha- 
los between the two that have similar masses (Evans et 
al. 2000) can be due to M31's higher initial densities by 
which star formation and chemical enrichment proceeded 
more efficiently before reionization compared with the 
Galaxy. We also suggest that the physical reason for 
the observed possible difference in the mean metallicities 
of stellar halos between M33 ([Fe/H] ~ —1.3; Brooks et 

al. 2004) and the LMC ([Fe/H] 1.5; Borissova et al. 

2004) might well be the same as above for the Galaxy- 
M31 difference. 

Although observations suggest that structural proper- 
ties and metallicity distributions of galactic stellar ha- 
los in the LG could be quite diverse (e.g., Durrell et al. 
2001; Alves 2004; Brooks et al. 2004), it remains to- 
tally unclear whether such diversity can be seen in other 
galaxies beyond the LG. Our simulations showed that 
physical properties of stellar halos in galaxies can be di- 
verse, partly because initial protogalactic conditions of 
a galaxy (e.g., the mass density) can determine to what 
extent reionization can influence star formation histories 
of subgalactic clumps that form the galaxy. In order to 
reveal radial structures of stellar halos in galaxies be- 
yond the LG and possible influences of reionization on 
the halo formation in these galaxies, we need to conduct 
wide-field imaging for global light distributions in galax- 
ies and measure the surface brightness distributions of 
the halos down to 28 — 30 mag arcsec~ 2 . 

For example, if we try to cover the 60kpcx60kpc halo 
region for a galaxy with the distance of ~ 10 Mpc, we 
need 20 x 20 wide-field imaging by a large telescope 
than can detect low-surface brightness stellar halos. Such 
wide-field imaging naturally enables us to reveal some 
stellar halo properties possibly influenced by reioniza- 
tion, such as the outer sharp edges in the surface bright- 
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ness distributions and the steep radial luminosity pro- 
files. Statistical studies of stellar halo properties by wide- 
field imaging for galaxies in different environments (e.g., 
clusters vs fields) can also help us to understand different 
degrees of reionization effects on the stellar halo forma- 
tion in different environments, which could result from 
possible inhomogeneity in reionization processes (e.g., 
Benson et al. 2001; Wyithe & Loeb 2004). This type 
of observations might well be formidable tasks for 8m- 
telescopes with wide-field imaging facilities: for instance, 
a few hour exposure with Subaru Suprime-Cam would 
reach as deep as 29 mag arcsec -2 in the £?-band with 
sufficient S/N ratios. Systematic observational studies 
of stellar halos by such telescopes will enable us to dis- 



cuss differences in reionization influences on stellar halo 
formation between different galaxies and thus to under- 
stand better the origin of galactic stellar halos. 
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Fig. 1. — Projected radial density profiles (E B (R)) of stellar halos for the standard models with z rc i on = 15 (open circles), z rc i on = 10 
(open square), and z re ion = 6 (open triangle). For clarity, T, B (R) normalized to the central value in each model is shown. For comparison, 
the power-law density profiles with S s (ij) cx R a for a = —1.0, —1.5, —2.0, and —2.5 are shown by dotted lines. Note that models with 
higher z re ; on show steeper profiles. The apparent bump at R ~ 6kpc for z rc i on = 15 is due to a stellar substructure. 
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Fig. 2. — The 2D distribution of surface brightness fiB (mag arcsec - 2 ) for the two extreme cases with z re i on =15 (left) and z r eion=6 
(right) for the standard models. This figure illustrates how z rc ; on affects the observed fig distribution of stellar halos. To derive these 
fig distributions, we divide the 60kpc X 60kpc halo region of a simulation into 50 X 50 meshes and thereby estimate fig for each mesh 
according to the photometric model described in the main text. We also smooth out the 2D distributions by using a Gaussian kernel with 
the smoothing length of 6 kpc. 
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Fig. 3. — The same as Figure 3 but for the two models with the same z re i on (=10) yet different initial over denisitcs: Si = 0.26 (left) 
and <5i = 0.52 (right). 



